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We study the unitary propagation of a two-particle one-dimensional Schrodinger 
equation by means of the Split-Step Fourier method, to study the coherent evo¬ 
lution of a spatially indirect exciton (IX) in semiconductor heterostructures. The 
mutual Coulomb interaction of the electron-hole pair and the electrostatic potentials 
generated by external gates and acting on the two particles separately are taken into 
account exactly in the two-particle dynamics. As relevant examples, step/downhill 
and barrier/well potential profiles are considered. The space- and time-dependent 
evolution during the scattering event as well as the asymptotic time behavior are 
analyzed. For typical parameters of GaAs-based devices the transmission or reflec¬ 
tion of the pair turns out to be a complex two-particle process, due to comparable 
and competing Coulomb, electrostatic and kinetic energy scales. Depending on the 
intensity and anisotropy of the scattering potentials, the quantum evolution may re¬ 
sult in excitation of the IX internal degrees of freedom, dissociation of the pair, or 
transmission in small periodic IX wavepackets due to dwelling of one particle in the 
barrier region. We discuss the occurrence of each process in the full parameter space 
of the scattering potentials and the relevance of our results for current excitronic 
technologies. 
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I. INTRODUCTION 


Excitonic states are the fundamental below-gap resonances in semiconductors heterostruc¬ 
tures, where stability and strength are strongly enhanced by quantum confinementj - 1 Spa¬ 
tially indirect excitons (IXs) are photo-excited correlated electron-hole pairs with the two 
charges localized in different segments of a nano-material. For example, IXs can be easily 
and selectively excited in semiconductor coupled quantum wells (CQWs), where an electric 
held applied perpendicular to the plane of the quantum wells separates positive and negative 
charges in different layers .®" 4 Since optical recombination depends on the electron-hole over¬ 
lap, this extends the intrinsic lifetime of excitons by orders of magnitude, from nanosecond^® 
to microseconds,^ 7 which allows for the formation of new types of quantum condensates.®® 

In CQWs with an average inter-well separation d, IXs are two-dimensional quasi-particles 
in, say, the (x, y) plane carrying a permanent electric dipole of order — ed along the growth 
direction z. Beyond any change in the binding energy implied by the particle separation, 
an electric held F z shifts the energy of IXs by U — —edF z . Therefore, although electri¬ 
cally neutral, IXs couple to any gradient of F z (x,y ) parallel to the planes .® 11 Such energy 
landscapes U(x, y) may be generated by imperfections and defects, but they can also be 
engineered by properly designed split gates, as in traditional nano-electronics, to modu¬ 
late the perpendicular held in the planes. Within their extended lifetime IX can be moved 
over several micrometers ,® 12 while recombination can be controlled by switching oh the per¬ 
pendicular heldf 13 broadening their application to an entirely new class of opto-electronic 
devices. In fact, in these devices data carried by photons can be stored in IXs, elaborated as 
in charge-based nano-electronics, and finally released back to photons by switching oh the 
perpendicular electric held at the desired time, all processes being integrated in the same 
semiconductor segment. Electronics based on excitons is called excitronics. Since gating 
can be modulated at the GHz frequency, that is orders of magnitude shorter than the IXs 
intrinsic lifetime, excitronic systems set themselves as a natural bridge between optical data 
communication and nano-electronic devices.®® 

Several excitronic functions have been demonstrated such as fast data storage,®® ac¬ 
celeration with electrostatic ramps 718 and interdigital devices ,® 9 held ehect transistors.® 
Most devices work in a dihusive regime, similar to traditional electronics. Therefore, IX are 
usually described theoretically as a classical gas of interacting dipoles .® 9 20 However, trap- 
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ping of single IXs has been recently demonstrated,^- 1 opening to the realization of quantum 
excitronic devices, similarly to single electron transistors (SET) in semiconductor quantum 

dotsPl 

Clearly, for perspective devices it is important to determine the behavior of IXs in different 
regimes, in order to optimize the performance. For example, large dipoles may improve 
acceleration and operation rate, but ensuing smaller binding energies may cause dissociation 
and loss of the stored information. IXs are indeed polarizable quasi-particles, with both 
center-of-mass and internal dynamics. Since in typical systems applied electrostatic fields 
in the CQW planes and the mutual Coulomb interaction of the pair are both in the few 
meV range, the quantum dynamics of IX under the influence of accelerating potentials is 
a genuinely two-body process, and the internal quantum dynamics cannot be neglected 
in general. Although a fundamental problem, the quantum description of such two-body 
scattering process has been little discussed in literature.^ 

In this paper we use a unitary propagation scheme to investigate theoretically the exact 
coherent dynamics of a single IX wavepacket describing scattering against different types of 
electrostatic potential landscapes, namely, steps/downhills and barriers/wells, using a min¬ 
imal ID model of a IX. In this sense, this is a two-particle generalization of the textbook 
description of one particle bouncing against potential steps and barriers. However, since the 
electron and the hole reside in different layers, in typical devices they are subject to differ¬ 
ent potentials. Therefore we have extensively investigated the parameter space separately 
for the two particles, analyzing the time-dependent dynamics during scattering as well as 
the long-time behavior of the scattered wavepacket. In addition to regions where IXs are 
reflected or transmitted almost as rigid objects, we identify experimentally relevant regimes 
where scattering is a genuine quantum two-body process. In these regimes, transmission 
or reflection of the wavepacket may be accompanied by excitation of the internal degrees 
of freedom, dissociation of the pair, or transmission in small periodic wavepackets due to 
dwelling of one particle in the barrier region. 

In the next section (Sec. 0 we set a ID model of a IX and we briefly outline the nu¬ 
merical method used for propagation. In Sec. |III| we analyze the space- and time-dependent 
dynamics of IX wavepackets scattering against potential steps/downhills and barriers/wells. 


In Sec. IV we discuss the result of scattering at asymptotic times in the full parameter space 
of scattering potentials. Our results are discussed in Sec. [V] in connection with experimen- 
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tally relevant regimes for current technologies. Details of the numerical method are given in 
Appendix [A} 

II. THEORETICAL APPROACH 

A. The Hamiltonian 

To study the dynamics of IXs, we use a ID model of an electron-hole pair propagating 
under the effect of the two-particle mutual Coulomb interaction and external electrostatic 
potentials coupling to the electron and hole separately. The two particles are tight to 
separated parallel channels at a distance d. Accordingly, we consider the following electron- 
hole Hamiltonian 


H = f + U c + U ext , 

(i) 



T - Pe + Ph , 

(2) 

2m e 

jj 1 e 2 /(47re 0 ) 

e r y/(£ e - X h ) 2 + d 2 ’ 

(3) 

U e xt Uefaeit) 

(4) 


Here p e (h)iX e (h), m e(h) are the ID linear momentum operator, the coordinate operator and 
the effective mass of the electron (hole), respectively. Uq is the electron-hole Coulomb 
interaction in a medium of relative dielectric constant e r . U e {h)(x e (h)->t) represents the one- 
body electrostatic potential which couples to the electron (hole) coordinate x e (h)- Note that 
in our model the internal dynamics of the electron-hole pair is described in full, i.e., we 
do not assume any rigid exciton approximationjZH2S nor the two particles are bound by 
construction: scattering against an external potentials may indeed result in excitation or 
dissociation of the pair. On the other hand, we neglect any lateral extension of the quantum 
states. Energy gaps induced by lateral confinement in a quantum well or wire (tens of 
meV) far exceed, and therefore decouple from, the small binding energy of IXs (few meV). 
Hence, the main effect of a lateral extension L is to regularize the Coulomb interaction at 
short distances to 1/L, 26 giving just a small renormalization of the parameter d. Finally, 
note that in U e (h)(x e (h),t) we indicated a possible explicit dependence on time t. Indeed, 
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the numerical methods which we discuss below also apply to the case of time-dependent 
potentials, although in this work we limit ourselves to stationary ones. 
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Figure 1. Schematics of the ID model assumed in the present work. 


Below we shall discuss simulations performed with two paradigmatic classes of electro¬ 
static potentials, namely, potential steps/downhills and potential barriers/wells. In SecjVjwe 
discuss the experimental relevance of these idealized potential profiles. A smooth potential 
step in the electron (hole) coordinate is described as a Fermi-like profile 


U e (h) (%e(h)) 


jj e i h )___ 

0 l _|_ e -(* e (h)-&)/“ ’ 


(5) 


x e (h) = b is the position of the step, U e ^) (b) = \U^ h \ For a particle hitting from the 
left, Uq^ > 0 represents a repulsive potential energy step, while U^ h ' 1 < 0 represents an 
attractive potential energy downhill, a is a smoothness parameter. Using a smooth potential 
avoids inaccuracies in numerical work, particularly when numerical derivatives need to be 
computed, and gives a better model of realistic potential modulations. 27 For simplicity, we 
assume common values a and b for electrons and holes. 

A potential barrier is described by a double Fermi-like profile 

U e(h) (x e{h) ) = U e 0 W x 


x 


1 

X _|_ e -(*e(Jv)-&l)/“ 
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1 

X -|_ e + he(t)-i>2)/a 



( 6 ) 


Here, b\, b 2 are the mid-potential positions for the rising and descending parts of the barrier, 
respectively. U [ ^ > 0 represents a potential barrier, while < 0 represents a potential 
well. 

In the absence of any external potential, namely for a free IX, it is natural to adopt the 
center of mass (CM) and relative coordinate system 


Y _ rn e x e +m h x h 
m B +m h 


X = X e ~X h 


(7) 
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Figure 2. Colormap of the potential landscape for an electron-hole pair in the presence of a barrier 
potential for the electron and a well for the hole with a = 2nm, b\ = 150 nm, 62 = 190 rim. Uq = 
+3.0 meV, t/g = —2.0meV as in Eq. §■ (a): (x e , Xh) representation, (b): (A, x) representation. 
Also indicated is the initial wavepacket, in arbitrary units, used in the simulations as described in 


Sec III The arbitrary origin is set at the center of the initial IX wavepacket. 


since the free Hamiltonian Ho = T + Uc splits into Hq = Hcm + H r where 

p2 


Hcm — 


2 M’ 


( 8 ) 


and 


1 eV( 4 «o) (9) 

2 m e r yjx 1 + d 2 

Here P, p are the CM and relative momenta, respectively, while M = m e + rrih and 
m = m e mh/M are the total and the reduced masses. The free IX wave function can be 
factorized as 


*(X,x,t) = $ K (X)Mx)e~ iEt/ \ (10) 

where the functions < I ) k(X") and <p n (x) are eigenfunctions of Hcm and H r with eigenenergies 
Ecm,k and E r)H , respectively. The CM component is labelled by the CM wavevector K and 
the relative component by the integer quantum number n. The zero of energy is set at the 
energy gap of the quantum well material, i.e., the energy for generating the two free particles 
at rest. Therefore E = Ecm,k + E r ^ n is the total energy of the IX. 

In the presence of an external potential there is no general solution in closed form, since 
coordinates (X, x) are coupled by U ex t- On the other hand, had we used the ( x e , Xh ) co¬ 
ordinates which are separable in U ex t, they would have been coupled by Uc- In principle, 
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the coherent propagation of a IX wavepacket can be carried out in both coordinate frames. 
However, as we aim at studying scattering events, we find the (X,x) representation more 
convenient to analyze the result at asymptotic times. Using the (X,x) set of coordinates 
turns out to be also beneficial from the numerical point of view. In fact, electrons and holes 
are usually in the same region of space, due to Coulomb attraction. As we use a real space 
representation, it is sufficient to describe small regions in r, that is r e ~ ry, rather then 
the full (x e ,Xh) space. 

Comparing the potential Uc + U e + Uh hi the two different coordinate systems may be 
useful to interpret the results of the following sections. In Fig. [2] we show the total potential 
in the case of a narrow barrier for the electron and a well for the hole, together with the 
square modulus of the initial IX wave function. The broad negative strip represents the 
Coulomb potential, while the narrow straight lines are the electron and hole single particle 
potentials, U e and Uh- In the (x e ,Xh) coordinate frame [Fig. |i](a)] U e and Uh are represented 
by straight stripes parallel to the axes, while the Coulomb interaction couples the two 
coordinates, being maximum along the diagonal x e = xy. On the contrary, in the (X,x) 
coordinate frame [Fig. |2j[b)] the Coulomb interaction only depends on the relative coordinate 
x, while the external potentials U e and Uh couple the (X, x) coordinates and therefore are 
represented by the two oblique stripes, the one with negative (positive) slope acting on the 
electron (hole) alone. Note the difference in the absolute values of these slopes, ensuing from 
the very different effective masses. 

B. Numerical propagation 

The quantum propagation of the IX is performed on a discrete homogeneous time grid, 
evolving the state between two consecutive times t and t + A t by applying the evolution 
operator U{t + A t , t ) to the wave function: 

T(x,t + A t ) = U{t + A t ,f)T(x,f). (11) 

The evolution operator is applied by using the Split-Step Fourier (SSF) method, which is 
numerically exact for vanishing A f . 28 29 Although this method is commonly exploited for 
the propagation of a single-particle wavepacket in a multi-dimensional physical domain, in 
Eq. we set x = (A", x). i.e, we interpret the two degrees of freedom of the system as the 
CM and relative coordinates of the electron-hole couple. 
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The SSF method relies on the Suzuki-Trotter expansion 30 ^ to factorize the evolution oper¬ 
ator as a product between exponential operators, each of which is diagonal either in direct or 
in Fourier space (details are given in Appendix [AJ. For this reason, this numerical method 
relies on a massive use of Fast Fourier Transformation (FFT), this resulting into a high 
efficiency and a relatively low computational cost in comparison with other methods. The 
use of FFT also implies periodic boundary conditions in x and X. However, the simulation 
domain is chosen sufficiently large that the wave function does not reach the boundaries in 
our analyses. Even if in this work we have been only concerned about time-independent 
potentials, the SSF method is also suitable for time-dependent ones. 


C. Initial state 

The time-dependent analysis starts from the choice of a realistic initial state. In a photo¬ 
generated exciton gas, the interaction with phonons and other scattering mechanisms allows 
for exciton relaxation to the ground state (j)o(x) within a few fs, i.e., on a much shorter 
timescale than photo-recombination lifetime. Therefore, at t = 0 we consider IXs to be in 
their ground state (fto(x) and we choose 


V{X,x,t = 0) = x{X)<f> 0 {x). 


( 12 ) 


where 0 n (a;), already defined in Sec. [TTj, are computed numerically via finite-difference diag- 
onalization, and where y (X) is a minimum uncertainty wavepacket 


x(x) 



' (X-X 0 )* - 

. . 


exp [iK 0 X], 


(13) 


centered in the initial CM position X 0 with width ax, and propagating with an average CM 
momentum 


Ko = 



Ecm, 


( 14 ) 


where Ecm is the most probable CM kinetic energy of the IX with total mass M. 








D. Analysis of the time-dependent dynamics 


In order to analyze quantitatively the time-dependent dynamics of the IX, we compute 
the quantum probability 

Pn(t)= If dxdX\q(X,x,t)\ 2 . (15) 


By choosing the integration domain, fl, as a suitable sub-domain of the simulation region 
(details are given in the following sections), we can quantify the transmission, reflection, 
and dissociation probabilities. 

Furthermore, we analyze excitations of internal degrees of freedom of the IX, which may 
take place due to the coupling between the internal degrees of freedom and the scattering 
potential, projecting the wavepacket on the eigenfunctions <p n {x) of the free Hamiltonian 
H 0 . Therefore, at any desired time t we comput^ 31 


c n (X,t)= / <p n (x)^(X,x,t)dx, 


(16) 


where Q is, as in Eq. (15), a suitable sub-domain. This quantity is further squared to 


eliminate complex phases and averaged over the CM coordinate, i.e. 


(|c,i(A",f)| 2 > x = J \c n (X,t)\ 2 dX. 
Moreover, we normalize the obtained coefficients, namely 

(|c„(V,«)| 2 ) x 


(17) 


k„r>(«) = 


Em(l C m(A',«)|%' 


(18) 


III. SPACE- AND TIME-DEPENDENT DYNAMICS 


We performed extensive numerical investigations of the wavepacket propagation in a 
broad range of the potentials intensities, Uq and Uq, considering both the attractive and the 
repulsive case. All simulations have been performed with typical GaAs parameters, as this is 
the material of choice for IX generation and propagation in CQW structures: m e = 0.067mo, 
rrih = 0.45 mo, where mo is the free electron mass, and e r = 12.9, corresponding to an effective 
Bohr radius as = 11.7 run . All simulations are performed with d — 17nm. 21 The eigenvalues 
of the IX free ID Hamiltonian H r obtained by straightforward numerical diagonalization are 
reported in Tab. [T| 
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n E rtU ( meV) 


0 

-4.63 

1 

-2.18 

2 

-1.26 

3 

-0.80 

4 

-0.55 

5 

-0.40 


Table I. Lowest six out of the eight negative eigenenergies of the free IX relative-coordinates Hamil¬ 
tonian Eq. ®> for the parameters used in the simulations (see text). 


The initial wavepacket has been chosen as described in Sec. |II C| with a most probable CM 
kinetic energy of Ecm = 0.5 meV. Typical grid point densities of 2.5 and 0.75 points/nm 
have been used for the x and X coordinate, respectively. Sharp potential profiles have been 
used with a = 2 nm. Propagation has been usually performed up to 40 picoseconds, with a 
time step A t = 20 fs, with (x,X) £ [— 0.6,0.6]/mi x [—1.5,1.7]/rm. 


A. Potential Step 

In this subsection we analyze the scattering against potential steps or downhills, as de¬ 
scribed by Eq. (|5]). It turns out that there are four typical phenomenologies, which we 
discuss below. 


1. Total reflection 

In Fig.g a) we show a typical evolution with Uq > Ecm and Uq = 0. This case resembles 
the textbook example of a particle reflected by a potential step, with interference fringes 
forming between incoming and reflected parts of the wavepacket. After the process is com¬ 
pleted, the square modulus of the wavepacket takes the original gaussian shape, spread in the 
CM coordinate according to the free single-particle propagation, with the standard deviation 
which can be obtained analytically to be crx(t) = crx(t = 0)y/l + {th/[2Mox(t = 0)]} 2 . 32 
This is the general behavior when Uq + Uq > Ecm■ Although it is in principle possible 
to transfer kinetic energy from the CM to the relative motion and excite IX internal lev- 
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els, for the present parameters this does not happen, since the CM kinetic energy Ecm is 
substantially smaller than the lowest excitation gap E r ^ — E r $ = 2.46 meV of the relative 
coordinate. 
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Figure 3. Evolution of a IX wavepacket scattering against a potential step/downhill set at b = 
150 nm. The origin of the CM coordinate is set at the center of the t = 0 wavepacket. In all cases, 
only one particle is accelerated: (a) Uq = 0.0 meV, Uq = +2.0 meV. (b) Uq = —4.0 meV, Uq = 0.0 
meV. (c) Uq = —5.0 meV, Uq = 0.0 meV. In (b) and (c) labels mark the integration domains used 


in Eqs. (16),(15) (see text), limited between dashed lines corresponding to x e = Xh = b (slanted 
lines) and x = ±15 as (vertical lines). A short animation of the continuous time evolution in (a) is 
included in the Supplemental Material. 33 
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2. Transmission 


In Fig. [3j^b) we show the complementary situation of an electron undergoing the accel¬ 
eration by a downhill, while the hole does not experience any potential. For this particular 
case, almost all the wavepacket is transmitted, as one would expect for a single particle 
scattering. However, here the transmission process is more complex. At intermediate times, 
for ps 12 ps, part of the wavepacket is strongly displaced from x — 0, that is the electron 
moves far from the hole, while the CM remains almost still (due to the larger hole mass, the 
CM position corresponds nearly to the hole coordinate). In this transient, the electron is 
accelerated, while the hole is not. At a later time, due to Coulomb interaction, the electron 
drags the hole along. In a classical picture, a wide oscillatory motion of the two particles is 
activated some picoseconds after the collision. In quantum terms, this corresponds to the 
excitation of higher quantum levels of the excitonic complex, which can be seen from the 
complex shape of the transmitted wavepacket in Fig. |3](b). 



Figure 4. Evolution 


of (|c n | 2 )(i) [see Eq. (18)], with H = A [see Fig. |3jb)]. 


To analyze quantitatively the evolution, in Fig. [4] we show the projection coefficients, 
(|c n | 2 ) (t). The average in the A" coordinate is performed over the transmission region 
A indicated in Fig. |3](b). The scattering starts at t ps 3 ps. Higher internal levels are 
increasingly excited and the scattering process is completed after ps 20 ps. The largest part 
of the wavepacket (ps 40%) is in the first excited (n = 1) state, while another ps 40% is 
almost equally contributed by states n = 0and2. All contributions sum up to ~ 96% of 
the wavepacket, which is the transmission probability. Note that different excited states are 
activated at different times and the most advanced part is of the n = 0 character. 
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3. Dissociation 


Similarly to excitation to higher IX states, dissociation into an unbound electron-hole 
pair is possible if 

(E r , o + E cm ) ~ (U e 0 + U%) > 0, (19) 

that is, an amount of energy exceeding the IX binding energy Eb = —E r 0 can be transferred 
to the relative degree of freedom. In our simulations, since Ecm Eb we can observe 
dissociation only with — (Uq + Uq) <4.1 meV. In our numerical evolution, dissociation 
occurs when the wavepacket asymptotically departs from x = 0 after scattering. In order to 
quantify this phenomenon, we compute Pq= a 3 (£), i.e., we integrate the probability density 
in a region far from x = 0 [see Fig. |3](c)] , and we define the dissociation probability 

Pdiss = JJ dxdX\V(X,x,t > too)! 2 . (20) 

where t^ is a sufficiently large time. 

In Fig. g c) a part of the wavepacket is pulled off to the right. Here Pdiss ~ 30%, as shown 
in Fig. §. Note that in this particular simulation Uq = 0, while Uq > Eb- Therefore, the 
electron acquires a high amount of energy from the downhill potential, this exciting the 
internal dynamics of the IX which eventually dissociates. This is similar (with role of the 
electron and the hole inverted) to Fig. |3](b) , where, however, Uq < Eb and the IX is excited 
but does not dissociate. 

4- Dwelling nearby X = b 

As a last interesting phenomenology arising from the two-body dynamics, we analyze 
cases in which part of the exciton wavepacket dwells in the proximity of X = b , i.e., near the 
step, for a long time interval. This occurs when both the electron and the hole experience a 
strong potential, but with opposite sign, a step for the electron and a downhill for the hole or 
vice versa. Figures |6](a),(b) shows the two cases. Dwelling may occur in region Ci [Fig. |6](a)] 
or C 2 [Fig. [6](b)] , depending on the electron/hole potentials. A semiclassical interpretation 
of the process is given in Fig. |6](c) with reference to the case shown in Fig. [6](b). As the IX 
wavepacket hits the scattering potential (1), the electron is accelerated along the downhill, 
while the hole does not overcome the repulsive step despite the Coulomb attraction from 
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Figure 5. Evolution of Ffo(i) [see Eq. (15)], with 0 coinciding with regions A 2 , A 3 , or the central 
region B indicated in Fig. [3][c). The probability over the region Ai is vanishing and it is not shown. 


the partner particle (2). Then, the electron, which is lighter, is pulled backwards by the 
Coulomb interaction (3). Part of its wave function is transmitted by the now-repulsive step, 
joins the hole at a small kinetic energy, and the whole IX is reflected backwards (4), while 
part of the electron wave function is reflected from the step, dragging the hole along, and 
the exciton as a whole is transmitted forward (5). However, there is a probability that, 
after reflection, the electrons does not pull the hole far from the barrier. In this case, the 
electron perform one more oscillation, and the process starts over again. Therefore, in such 
a process reflection and transmission occur as a periodic pumping of smaller wavepackets. 
The description of the process in Fig. [6|a) is analogous with the roles of the electron and 
the hole inverted. 

B. Potential Barrier 

We now investigate the scattering from a finite range potential, Eq. ©• A barrier/well 
for each carrier is present in a region 6 2 — &i = 40 nm wide. Note that, while the width 
of the barrier/well has an obvious effect on the scattering event, its absolute position is 
unimportant, as far as the two parameters b\ and b -2 are the same for the two particles. In 
the present case quantum tunneling becomes possible, as in usual single-particle scattering. 
On the other hand, in the asymptotic regions Uq = Uq = 0. Therefore, it is not possible 
to excite internal levels or to have free particles, since, as discussed in the previous section, 
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Figure 6. As in Fig.[3]with (a) Uf = +3.0 meV, Uf} = —2.0 meV.(b) Uq = —3.5 meV, Ujf = 1.5 meV. 
(c) Semiclassical interpretation of the scattering process shown in (b). Dashed lines represent 
equations x e = b and Xh = b. 


the CM kinetic energy which could be transferred to the internal degrees of freedom is much 
lower than the excitation gaps of the IX and, a fortiori, its binding energy. However, an 
oscillatory motion or even the dissociation of the IX may occur if one of the particle is 
captured by its well, transferring a large energy to the partner particle. 

In Fig. [Tj^a) a wavepacket hits a high barrier for electrons and a weaker well for the 
hole and tunnels with « 19% probability. Note that the hole dwells into its potential well 
in the initial stages of the process. Dwelling may also lead to periodic emission of the 
wave function, as for the step/downhill case, which is shown in Fig. [T^b) and sketched in 
Fig0 c). When the IX hits the barrier/well potential (1), the electron may be captured by 
its well, while the hole starts to oscillate around the well position (2). At each oscillation, 
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the wavepacket located in the well region is partially transmitted and partially reflected as 
a IX bound state (3). These oscillations are highlighted by the step-like evolution of the 
transmission coefficient, shown in Fig. [8j 
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Figure 7. Scattering through potentials barriers/wells with b\ = 150 nm, 62 = 190 nm. (a) Uq = 
3.5 meV, Uq = —1.5 meV. (b) Uq = —4.0 meV, Uq = 0.0 meV. (c) Sketch of the periodic emission 
of the IX wavepacket. Short animations of the continuous time evolution in (a) and (b) are included 
in the Supplemental Material . 33 


IV. SCATTERING PHASE DIAGRAM 

The space- and time-dependent evolutions discussed in the previous section are selected 
examples in a large set of simulations, to highlight how the scattering process takes place 
in different regimes. In this section we summarize the outcome of the scattering events at 
asymptotic times in the full parameter space Uq and Uq. 
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Figure 8. Transmission coefficient calculated from Eq. (15) with Q = A [see Fig. [7](b)]. 


For each simulation we computed transmission, reflection, and excitation probabilities 
at t = 40 psp3 that is, at a sufficiently large time for the scattering process to be com¬ 
pleted. These coefficients are used to identify distinct phenomenologies taking place at 
different values of the potential parameters Uq and Uq. For each class of potentials, either 
steps/downhills or barriers/wells, we summarize our results in a phase diagram as a function 
of the potential total intensity (Uq + Uq) and potential asymmetry (Uq — Uq) between the 
two particles, rather then Uq and Uq separately, as the former quantities are those driving 
the distinct behaviors. A common energy scale for all simulations is set by the kinetic energy 
Ecm — 0.5 meV. In the phase diagrams to be discussed below, it is convenient to report 
results in terms of the dimensionless ratios (Uq + Uq)/Eqm, (Uq — Uq)/Eqm ■ Note that 
for a rigid exciton, i.e., an exciton with frozen internal degrees of freedom, 24 -'the scattering 
outcome (i.e., transmission probability) scales with these ratios, similarly to the textbook 
single particle scattering. However, in the present case of a quantum system with an inter¬ 
nal degree of freedom, the latter couples to the CM kinetics and is affected by the external 
potentials. Therefore, here the scaling is only approximate: although still a convenient rep¬ 
resentation, one should keep in mind that different behaviors might be expected in a given 
region of the phase diagram if Ecu varied. This is particularly true for what concerns the 
phenomenologies arising from the internal dynamics, i.e. dissociation, excitation to higher 
levels, pausing/pumping nearby the external potential region. 
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(a) Steps/Downhills 



(U 0 e + U 0 h )/E CM 


(b) Barriers/Wells 



(U 0 e + U 0 h )/E CM 


Figure 9. Phase diagram for (a) steps/downhills and (b) barriers/wells simulations. The transmis¬ 
sion probability is displayed in grayscale. Simulations are performed every 1 meV in (Uq + Uq) 
and (Uq — Uq). The main phenomenologies are highlighted by superimposed letters and solid lines 
to mark approximate boundaries of different regions. Regions D: dissociation probability greater 
than about 5%. Regions E: excitation to higher internal levels greater than about 10%. Regions 
P: periodic transmission/reflection of the wave function due to dwelling. The positive (negative) 
slope dashed line indicates Uq^ = 0. 

A. Potential steps/downhills 

1. Transmission/reflection 

In Fig. § a), the transmission probability is reported as a gray scale colormap. Simulations 
are performed on a coarse grid with steps of 1 meV in (Uq + Uq) and (Uq — Uq), i.e., in steps 
of 2 units in the ratios (Uq ± Uq) / Eqm- This choice of the simulation space is a good trade¬ 
off between a sufficiently detailed insight and the computational cost. In Fig. |9](a) , each 
of the 242 coloured squares corresponds to a specific simulation running for about 8 hours 
(including postprocessing) on a high-end workstation. When the IX wave function leaves the 
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central region of the domain, i.e. the active region where the external potential may cause 
reflections or trigger excited internal states, the transmission probability is computed by 
means of Eq. (15). The integration region Q is indicated as A 2 in Fig. [3j^c) and corresponds 
to the probability of finding both the electron and the hole beyond the step (x e > b and 
Xh > b ) and less than 15a# apart (x < 15ob). 

As might be expected, transmission is large on the left part of the diagram, corresponding 
to electron and hole downhills, and drops on the right hand side, which correspond to 
repulsive steps for both particles, vanishing at (Uq + Uq)/Eqm — 1- Note that if the 
potential is symmetric for the two particles, i.e., (Uq — Uq) — 0, the transmission rapidly 
drops from 1 at (Uq + Uq)/Ecm — 0 to zero. However, for strongly asymmetric potentials 
(Uq — Uq)/Ecm\ > 5, that is in the upper and lower parts of the diagram, the transmission 
start to decrease at negative values of the total potential Uq + Uq and evolves smoothly to 
zero. This is because in these regions, one particle experiences a strongly attractive downhill 
potential, while the other is subject to a weaker but repulsive potential. Note also that the 
diagram is slightly asymmetric in the upper and lower parts. This is due to the different 
masses of electrons and holes, so that similar potential intensities, but with the role of the 
particles inverted, have a different effect on the transmission. 


2. Excitation/Dissociation 

Excitations of the relative motion eigenstates 4> n (x) are favored when both (Uq + Uq) < 0 
and (Uq — Uq) < 0 [lower left corner of Fig. |9^a)], which corresponds to a large downhill 
potential for the electron and a weaker potential step for the hole. Clearly, the situation 
may be reversed, with electron and hole potentials inverted, (Uq — Uq) > 0 (top left corner). 
However, in this case the excitation is typically weaker, so that regions of excitations and 
dissociation move to stronger (positive) anisotropies in Fig. |9](a). 

This asymmetry can be understood from inertial considerations: since the hole is one 
order of magnitude heavier than the electron, a step/downhill for the hole almost coincides 
with the potential acting on the CM coordinate [a different way to see this is to note that 
the potential step for the hole is nearly independent from the relative coordinate, compare 
Figsj3](a) with Figs|3](b),(c)], and the IX accelerates as a single, rigid particle. On the other 
hand, when an electron experiences an equivalent downhill potential, a larger part of the 
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external energy can be transferred to the internal relative dynamics. Dissociation can also 
take place, when the electron downhill is particularly strong. 35 


B. Potential barriers/wells 

Contrary to the step/downhill case, with a barrier/well potential it is not possible to 
have excitation of internal levels nor dissociation in the transmitted part of the wavepacket. 
As discussed in the previous section, however, for (Uq + Uq)/Ecm •C 0 and \U^ h) \ 
Eqm dissociation in the barrier region takes place if one particle (the electron or the hole, 
depending whether Uq or Uq is negative and intense) is captured into its well, while the 
other inertially continues its motion. 

Due to quantum tunneling, the transmission probability is finite also for (Uq+Uq )/Eqm > 
1. Note that tunneling is substantially anisotropic between the upper part (Uq > 0 and 
Uq < 0) and the lower part (Uq < 0 and Uq > 0) of the phase diagram. Again, this 
is clearly due to the difference between electron and hole effective masses, as for large 
anisotropies IX does not tunnel as a rigid complex of mass M, but rather through a complex 
two-body process. On the other hand, as one may expect from one-body tunneling, resonant 
transmission may take place. This is shown in Fig. [9](b) by the oscillatory behavior of the 
transmission coefficient, particularly in the upper part of the phase diagram. 

Finally, note that, at difference with the steps/downhills case where dissociation and 
excitation take place more easily for negative anisotropies (lower left part of the phase 
diagram) than for positive anisotropies (upper left part), here the opposite is true, since 
these phenomena happen by the trapping mechanism, which is much more favorable for 
holes. 

V. CONCLUSIONS 

We have discussed the unitary space- and time-dependent quantum dynamics of a single 
IX in a CQW system scattering against simple potential profiles. Our analysis suggests, in 
particular, that exciton scattering in typical devices is a genuine two-body process, substan¬ 
tially different and more complex than the single-particle scattering of a rigid exciton. Our 
results are based on a minimal ID model. Although strictly speaking this analysis applies to 
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ID channels (e.g., electrostatically defined coupled quantum wires), it allows to capture the 
main phenomenologies to be expected also in planar CQWs. Next we discuss the connection 
to realistic systems and perspective experiments. 

For the prototypical potential profiles considered here (but in a large range of parameters) 
IX scattering times turn out to be in the order of tens of ps. This is smaller than the 
lifetime of IXs in CQWs and also their coherence time in quantum wells, 36 and supports the 
possibility to engineer single-IX phase coherent excitronic devices. 

Having used typical material parameters and potential landscapes for GaAs-based struc¬ 
tures, we expect that the complex phenomenologies found in our simulations could take place 
in excitronic devices and should be taken into account in device engineering. Moreover, the 
complex time-dependent IX dynamics exposed by our results could be directly probed, for 
example, by collecting time-resolved micro-photoluminescence maps' 12 during the scattering 
process, exploiting the fact that optical recombination of IXs can be triggered at a desired 
time by switching off the vertical field F z ) 17 1 37 Luminescence maps could be calculatecP®® 
from |T(A", x = 0)| 2 , as provided, e.g., by our simulative method. Although we did not 
attempt an explicit calculation of optical signatures to be directly compared to experiments, 
since our model is strictly ID, for purpose of illustration we provide in the Supplemental 
Material 33 Uhe real-space and K-space evolution of T(X, x = 0). Note that in this respect 
single-IX devices could offer a unique possibility to probe the quantum two-body scattering 
dynamics by optical means. In realistic 2D devices, scattering times could be larger than 
those found in our simulations, due to the larger phase space. Time-resolution can thus be 
less demanding. 

Potential profiles as those employed in the present study are usually induced by electro¬ 
static gates. On the one hand, typical top gates induce an electrostatic potential in the 
underlying 2D system with opposite sign for the two particles. Due to layer separation d, 
however, the potential generated by a top gate is different in magnitude for the two layers. 
For example, for a device of total thickness 1 /irn with 1 V applied to top gates the resulting 
voltage difference between the two CQWs 10 nm apart is of the order of 10 mV. Therefore, 
a common situation in excitronic devices is with large potential anisotropies, corresponding 
to the upper and lower sectors between the bisecting lines in Fig. [9j 

On the other hand, it is also possible to gate separately the two layers of the heterostruc¬ 
ture and engineer potential landscapes of the same sign and similar intensities for the two 
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particle^- 41 ! Potential profiles of the same sign for the two particles with few meV anisotropies 
correspond to left and right regions between the bisecting lines in Figs. [9] Scattering poten¬ 
tials of the same sign may also arise from monolayer fluctuations of the QW thickness which 
induce modulations of the lateral confinement energy in the meV range i.e, a steps/downhills 
potential profiles along the planes of the QWs. 
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Appendix A: The two-particle Split-Step method 


In order to study the time-dependent scattering dynamics of IXs, we use a numerically 
exact solution of the Time Dependent Schrodinger Equation (TDSE) with the Hamiltonian 
given by Eq. ([Tj) based on the Split Step Fourier method. 

The time propagation of a generic quantum state T(x, t), between two close time instants, 
t and t + A t , can be expressed by means of the evolution operator, U{t + A t ,t), as: 


T(x,t +A t ) = U(t + A t ,t)T(x,t), 


where: 


+ At, t ) — 

= Texpj-^-^ (f + U( T )) dr\ « 

K Texp{-i(fA t + C /( t + ^)A,)}. 


(Al) 


(A2) 


Here Texp{} is the so-called time ordered exponential operator. In the last line the ‘rectangle 
integration method’ has been used for integrating the potential energy. The last line in 


Eq. (A2) is exact for time independent potentials. 


22 



By using the Suzuki-Trotter formula 30 in order to expand the exponential of the sum of 


operators in Eq. (A2) as a product of exponential operators, we obtain 

U(t + A t, t ) = 

e -i^+^)# e -|TA te -|C7( t+ ^)^ + 0(A 3 } _ 


(A3) 


We now exploit the fact that the operator T is diagonal in Fourier space, while the 
operator U is diagonal in direct space. The basic idea behind the Split Step Fourier (SSF) 


method is to exploit the evolution operator factorization, Eq. (A3), in order to act on the 


state function in the specific representation in which the operator is diagonal. At each time 
step the SSF algorithm performs the following operations: 

1. Multiplication of the wave function in coordinates representation by the rightmost 
potential term, for half time step: 


T(x, t) —)• T'(x) = e 




T(x,t). 


(A4) 


2. Computation of the Fourier transformation of the state T'(x): 

T'(p) = [ dxe- ip x T'(x). 


(AS) 


3. Multiplication of the wave function in momentum space by the kinetic energy term 
for a full time step, A t : 


*'(P) —► *"(p) = e 




*'(P). 


(A6) 


4. Calculation of the inverse Fourier transform of the wave function: 

1 


x = 


(2ny 


dpe +ip - x T"(p). 


(AT) 


5. Multiplication of the wave function by the leftmost potential term for another half 
time step: 

"*)3Vfx). (A8) 


T"(x) —> T(x,t + A t ) = 


The SSF is a unitary method which preserves the norm of the wave function. 28 We underline 


that in Eq. (A3) operators T and U are represented in their respective diagonal forms using 


the appropriate basis. Therefore, even though a spatial discretization of the domain (as 


23 






described in Sec. Ill) is used in order to represent the wavepackets and the potential U, 
no finite difference approach is needed to represent the kinetic operator, which instead is 
represented in the reciprocal grid implied by the discrete Fourier transform algorithm. 

The SSF method is commonly used to propagate a wavepacket representing one particle 
in a n-dimensional space, e.g. with x = (x, y) for n = 2. However, since it is an exact 
solution of the Schrodinger equation, there is no limitation as to the interpretation of the 
set of coordinates x. Here, we implement the method by propagating the two-body wave 
function in the x = (X,x) coordinates. Note that the momentum vector p = (P, p) has also 
components corresponding to the CM and relative momenta, respectively. 

Unlike the isotropic mass one-body case, the two coordinates are governed by different 
terms in the Hamiltonian, with very different effective masses, resulting in different frequency 
ranges for P and p. While typical energies for the relative motion are in the few meV range, 
the typical kinetic energy for the CM is in the fraction of an meV. It is this decoupling which 
makes particularly convenient to work in the CM and relative coordinates. 


Appendix B: Wavepacket in K space 


In Fig. 10 we show the A'-space representation of H/(A, x = 0) for the three cases (and 
the same selected times) presented in Figs. 3(a), 7(a) and 7(b). This quantity is directly 
related to the luminescence intensity 38 40 which is oc |\k(A, x = 0)| 2 . The highlighted region 
around the origin satisfies 


77 

K < -E gap , (Bl) 

which represents the A'-space radiative cone.® Here, n = sJEr is the refractive index of GaAs, 
c the speed of light, h the reduced Planck constant and E gap = 1.424 eV the GaAs energy 
gap. 

The videos movieK_f ig3a.mpg, movieK_f ig7a.mpg and movieK_f ig7b.mpg of the Sup¬ 
plemental Material 3 ® show the time evolution of |'F(X, 0)| 2 (top panel) and its A'-space 
transform (bottom panel) for the three cases of Supplementary Fig. 1. The K components 
that comply with Eq. (1) are also highlighted. The evolution shows that such optically-active 
components are substantially modified during the scattering of the IX. 

Moreover, the videos movieXx_f ig3a.mpg, movieXx_f ig7a.mpg and movieXx_f ig7b.mpg 
of the Supplemental Material 3 ®, show the continuous-time evolution of |T(A", x)\ 2 for the 
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Figure 10. A"-space representation of the wave function at four different times for the same IX 
scatterings shown in Figs. [3^a), [7](a) and[7]^b), in columns (a), (b) and (c), respectively. 

three scattering events also shown in Figs. 0a). 0a) and (7^b). 
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